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We analyse the effect of bounding the occupation number of bosonic field modes on the correlations 
among all the different spatial-temporal regions in a setting in which we have a space-time with 
a horizon along with an inertial observer. We show that the entanglement between A (inertial 
observer) and R (uniformly accelerated observer) depends on the bound TV, contrary to the fermionic 
case. Whether or not decoherence increases with TV depends on the value of the acceleration a. 
Concerning the bipartition AR (Alice with an observer in Rindler's region IV), we show that no 
entanglement is created whatever the value of TV and a. Furthermore, AR entanglement is very 
quickly lost for finite TV and for TV — > oo. We will study in detail the mutual information conservation 
law found for bosons and fermions. By means of the boundary effects associated to TV finiteness, 
we will show that for bosons this law stems from classical correlations while for fermions it has 
a quantum origin. Finally, we will present the strong TV dependence of the entanglement in RR 
bipartition and compare the fermionic cases with their finite TV bosonic analogs. We will also show 
the anti-intuitive dependence of this entanglement on statistics since more entanglement is created 
for bosons than for their fermion counterparts. 
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I. INTRODUCTION 

To analyse quantum correlations in non-inertial set- 
tings it is necessary to combine knowledge from differ- 
ent branches of Physics; quantum field theory in curved 
space-times and quantum information theory. This com- 
bination of disciplines became known as relativistic quan- 
tum information, which is developing at an accelerated 
pace . It provides novel tools for the analysis of the 

Unruh effect [13 - [2l| allowing us to study the behaviour 
of the correlations shared between non-inertial observers. 

Let us consider a bipartite system (Alice- Rob), in 
which one of the partners (Rob) is undergoing a uniform 
acceleration and therefore describing the world which 
Rindler coordinates. As pointed out in 0, [23| there 
are 3 possible bipartitions that can be considered when 
analysing entanglement in the previous setting; 1) The 
entanglement of the inertial observer with field modes in 
Rindler's region / (Alice- Rob AR), 2) The entanglement 
of the inertial observer with field modes in Rindler's re- 
gion IV (Alice-AntiRob AR) and 3) The entanglement 
between modes in regions I and IV of the Rindler space- 
time (Rob-AntiRob RR). Partitions AR and AR are 
especially important as these are the partitions in which 
classical communication is allowed (we will refer to them 
as CCA bipartitions from now on). 

In previous works some aspects about correlations 
comportment in the presence of horizons were disclosed. 
These works accounted for the radical differences be- 
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tween fermionic and bosonic entanglement behaviour in 
the presence of Rindler and event horizons, showing that 
the real cause of these differences is fermionic/bosonic 
statistics. This contradicts the naive argument that the 
differences come from the finite dimensional nature of the 
fermionic Hilbert space for each frequency as opposed 
to the built-in infinite dimension of the Fock space for 
bosons. 

However, in previous works it was pointed out that 
there exist effects of the dimension of the Fock space 
on the correlations in the bipartition RR. As it was 
shown in [23[ these effects clearly distinguish the case 
of fermions with a different Fock space dimension (Dirac 
vs Grassman scalar). It was an open question, however, 
to account for the precise effect of the dimension of the 
Fock space for each mode in the correlations behaviour 
of bosonic fields. In this work we will see that, through 
these effects, we can reach a deeper understanding about 
the behaviour of correlations in the presence of horizons. 

Nevertheless, in earlier studies it was pointed out that 
there exist effects of the dimension of the Fock space 
on the correlations in the bipartition RR. As it was 
shown in [23| these effects clearly distinguish the case 
of fermions with a different Fock space dimension (Dirac 
vs Grassman scalar). However, It is still an open question 
to account for the precise effect of the dimension of the 
Fock space for each mode on the correlations behaviour 
of bosonic fields. In this work we will see that analysing 
these effects we can reach a deeper understanding about 
the comportment of correlations in the presence of hori- 
zons. 

In an attempt to account for the dimensional effects 
mentioned above, we are going to study a finite dimen- 
sional analog to bosonic fields with a limited dimension 
of each frequency Hilbert space. This means engineering 
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a method to impose a maximum occupation number N 
in each scalar field frequency mode. The construction of 
a finite dimensional scalar field state for a non-inertial 
observer can in itself be problematic, thus an issue which 
will need to be tackled in order to conduct the proposed 
analysis. 

We will present results about entanglement of the CCA 
bipartitions that will strengthen the argument discussed 
in previous works (THHHIill on the capital importance of 
statistics in the phenomenon of Unruh decoherence and 
its role in the black hole information paradox. Specif- 
ically we will prove that the behaviour is fundamen- 
tally independent of the Fock space dimension. However, 
bosonic entanglement for AR is slightly sensitive to Fock 
space dimension variation, in opposition to what happens 
with fermions. 

We shall point out that those variations are of most 
interest themselves. Namely they again strongly oppose 
what is said in previous literature that Unruh decoher- 
ence degrades the entanglement quicker as the dimen- 
sion of the Hilbert space is higher. Instead, we will show 
that quantum correlations can be more or less quickly de- 
graded for different dimensions depending on the value 
of the acceleration. This completely banishes the former 
argument. 

We shall point out that those variations strongly op- 
pose once again what is said in previous literature in 
which it is argued that the Unruh decoherence degrades 
the entanglement quicker as the dimension of the Hilbert 
space is higher. Instead, we will show that quantum cor- 
relations can be more or less quickly degraded for differ- 
ent dimensions depending on the value of the accelera- 
tion. This completely banishes the former argument. 

We will also show how classical correlations between 
AR and AR are affected by the bound on the occupa- 
tion number. Specifically we will show that the effect of 
imposing a finite dimensional Fock space affects the con- 
servation law for mutual information found in previous 
works. We will compare this with the fermionic cases and 
will prove some results about mutual information to be 
completely universal. 

Furthermore, we will show remarkable results concern- 
ing correlations between modes in Rindler regions I and 
IV. We will see how they are ruled by bothstatis- 
tics and Hilbert space dimension. There arc differences 
and similarities between fermions and bosons concerning 
correlations RR. We will analyse the different bosonic 
cases comparing them with their Fock space dimension 
fermionic analogs, in order to comprehend the relative 
importance of dimensionality and statistics in the be- 
haviour of such correlations. 

This work is structured as follows. Section ITll provides 
a brief introduction about Rindler space-time and the 
perspective of an accelerated observer. In section IIIII 
we will construct bosonic bipartite entangled states with 
bounded occupation number. In section HVl we will anal- 
yse entanglement (subsection IIV A| and mutual informa- 
tion (subsection IIV Bp for all the possible bipartitions of 
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Figure 1: Rindler space-time diagram: lines of constant posi- 
tion z = const, are hyperbolas and all the curves of constant 
proper time t for the accelerated observer are straight lines 
that come from the origin. An uniformly accelerated observer 
Rob travels along a hyperbola constrained to region I 

the system and different values of the occupation num- 
ber bound N, comparing the results of finite dimensional 
bosons with the fermionic analogs. Finally, we present 
our conclusions in section IVl 

II. SCALAR FIELDS AS SEEN FROM AN 
ACCELERATED OBSERVER 

A uniformly accelerated observer describes the world 
by means of a set of Rindler coordinates [24j . We need, 
however, two different set of coordinates in order to cover 
the whole Minkowski space-time. These sets of coordi- 
nates define two causally disconnected regions in Rindler 
space-time. If we consider that the uniform acceleration 
a lies on the z axis, the new Rindler coordinates (t, x, y, z) 
as a function of Minkowski coordinates (t, x, y, z) are 

at = e az sinh(at), az — e az cosh(ai), x — x, y = y (1) 

for region I, and 

at = —e az sinh(at), az — —e az cosh(at), x = x, y = y 

(2) 

for region IV. As we can see from fig. I, although we 
have covered the whole Minkowski space-time with these 
sets of coordinates, there are two more regions labeled II 
and III. To map them we would need to trade cosh by 
sinh in equations (H},©. In these regions, t is a space- 
like coordinate and z is a timelike coordinate. However, 
considering such regions is not required to describe fields 
from an accelerated observer perspective P.l5lll5l.l23.l25j. 

The Rindler coordinates z,t go from -oo to oo inde- 
pendently in regions I and IV. It means that each re- 
gion admits a separate quantisation procedure with their 
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corresponding positive and negative energy solutions of 
Klein- Gordon equations , ipl~ and ipf. , ^k^- P° s ~ 
itive and negative energy solutions will be classified with 
respect to the future-directed timelike Killing vector in 
each region. In region I the future-directed Killing vector 
is 



The vacuum state for a fc-momentum mode of a scalar 
field seen from the perspective of an accelerated observer 
is 



|0) M = 



1 



coshi 



;tanh n r \n) 1 \n) IV 



(7) 



n=0 



dI t=Ft di+ Tt di = ai ~ Zdi 



■ids), 



(3) 



whereas in region IV the future-directed Killing vector is 

or = -di 

This means that solutions in region I, having time de- 



Notice that here and from now on we are dropping the k 
label as we are working under SMA. 

The Minkowskian one particle state results from ap- 
plying the creation operator to the vacuum state. Its 
translation to the Rindler basis is 



pendence ip k + ~ e with > 0, represent positive _ 1 

energy solutions, whereas solutions in region IV, having M cosh 2 r _ Q 

t.imp Hpnpnrlp-nr'P lhl^~ ^ p~ ikot w itV, U n ^> fl n r p flptnallv 



oo 

^ tanh' 



time dependence tpk + ~ e~ lkot with k > 0, are actually 
negative energy solutions since 9/ points to the oppo- 
site direction of d^. As I and IV are causally disconnected 
tp^y and ipl^ only have support in their own regions, 
vanishing outside them. 

Let us denote a T k , at k as the particle annihilation and 



i n rVn + l \n+l)j \n) IV 



creation operators in region / and a 



IV,ki a IV,k 



as the par- 



ticle/antiparticle operators in region IV. 

The bosonic operators satisfy the commutation rela- 
tions [a R k, a li' fc'] = ^RR'^fcfc'- The subscript R notates 
the Rindler region of the operator R = {/, IV}. All other 
commutators are zero. 

We can relate Minkowski and Rindler creation and an- 
nihilation operators by taking appropriate inner prod- 
ucts and computing the so-called Bogoliubov coefficients 
[!, HI HI, HI ■ This allows us to translate Minkowskian 
states into Rindler coordinates as we will present in the 
next section. For a scalar field, the Bogoliubov relation- 
ships for the annihilation operator of modes with positive 
frequency are 



where 



O'M.k — cosh r ai t k — sinh.ra\ v _ k 



tanh r 



(4) 



(5) 



III. LIMITING THE OCCUPATION NUMBER 

The Minkowski vacuum state for the scalar field is de- 
fined by the tensor product of each frequency mode vac- 
uum 



10) 



|o A 



(6) 



such that it is annihilated by a k for all values of k. 

Through this work we will operate under the single 
mode approximation (SMA) 0, In any case in [22[ 
we discussed the implications of such approximation and, 
for the purposes of this work, going beyond SMA will not 
give relevant differences. 



Now we will consider the following maximally entangled 
state in the Minkowsky basis 

|*) = i=(|0) M |0) M + |l) M |l> M ) (9) 

This is a qubit state which is a superposition of the bi- 
partite vacuum and the bipartite one particle state. 

For our purposes we need to limit the dimension of 
the Hilbert space. To do so we are going to limit the 
maximum occupation number for the Rindler modes up 
to N. 

As a first approach, one could consider bosons with a 
limited dimension of the Hilbert space by imposing oc- 
cupation number bounds on the states ([7]), ([5]) to model 
impenetrable bosons as done to illustrate some qualita- 
tive arguments in [23j . 

However, if we want to go beyond the qualitative ef- 
fects and do a completely rigorous analysis we would 
run into important problems. Namely, we would not be 
able to normalise both states ([7]), © simultaneously. In 
other words, the translation into Rindler coordinates of 
the Minkowskian creation operator applied on the vac- 
uum state do not preserve normalisation and it is no 
longer true that applying the annihilation operator to 
the one particle state we recover the vacuum of the the- 
ory. Therefore the second quantisation rules of bosonic 
fields would be ill-defined (for instance, problems would 
appear when applying the commutator to the one par- 
ticle state). As statistics is fundamental to explain the 
Unruh decoherence mechanism, a rigorous analysis would 
require that we consider the vacuum of our theory as ex- 
pressed in (J7J with an unbounded occupation number. 

Alternately, we will define finite dimension analogs to 
the vacuum and one particle states 



I°Jv)m 



1 



N 



cosh i 



tanh" r \n) j 



/IV 



(10) 



n=0 



N-l 



|!iv) M = 72~ / ] tanh'V \/n+ 1 \n+ l)j \n) IV 

cosh r n 

(11) 
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in which we have cut off the higher occupation numbers 
and thus these two states are not exactly the vacuum of 
our theory and the first excitation. Instead, they could be 
understood as approximations in which Rob is not able 
to notice occupation numbers larger than N. This sim- 
ple construct allows us to consider a bounded occupation 
number along with bosonic statistics. Therefore we can 
now disentangle the statistical effects from the ones de- 
rived from the dimensionality of the Hilbert space. This 
limiting our states to a subset of relevant modes is a 
similar procedure to the one used to carry out the SMA 

We will then consider the following entangled state in 
Minkowsky coordinates 



I*) 



1 



Cjv(r) 



(|0) M \°n)m + |1)m IMm) (12) 



in which the one particle state and the vacuum for Rob 
are substituted by the bounded occupation number ap- 
proximations. 

Notice that a factor 1/Cjv(r) must now be included 
as our occupation number cutoff implies that \0n) M and 
|1jv) m are not normalised. Its value is 

Cn{t) = ^{Qn\0n) m + {1n\In) m (13) 




2 - tanh 2W r ( tanh 2 r + 1 



N 



cosh 2 r 



(14) 



In the limit N —> oo, CV(r) -+ V% recovering the stan- 
dard scalar maximally entangled state ([§]). 

This normalisation will also guarantee that Unruh 
decoherence will not affect higher occupation number 
modes, in similar way as SMA restricts the decoherence 
to a single frequency. 

The density matrix for the whole tripartite state, 
which includes modes in both sides of the horizon along 
with Minkowskian modes, is built from (fT2|) changing to 
Rindlcr coordinates for Rob 



The bipartite systems are characterized by the follow- 
ing density matrices 



, n=0 



cosh 2 : 



\0n)(0n\ + ^-^(\0n)( 
cosh r \ 

)7T f —I— 1 
+ — n-l ln + 1 )( ln + 1 l 
cosh" r 



ln + 1 



tanh 2Ar r 
cosh 2 r 



\0N)(0N\ 



1 



Cjv(r) 



(22) 



'^tanh^r 



E 

n=0 



cosh 2 r 



| On) (Ora| 



Vn+ 1 



cosh 

71+1 



1 / 

— tanhrl |0 ra + 1) 
r V 



x (lra| + |lra)(0ra+ 1|) + - |lra)(lra| 



cosh r 
1 



^^|OA0(OA|}--^ 

cosh r C N (ry 



(23) 



p RR = 



N 

E 

m=0 



tanh n+m r , . , v-^ 1 tanh n+m r 



cosh 2 r 



|rara) (mm\ + 



n=0 
m=0 



x \Jn + lvm + 1 \n + 1 ra) (rra + 1 m\ 

where the bases are respectively 

\nm) = \n A ) M \m R ) I 
\nm) = \n A ) M \m R ) IV 
\nm) = \n R )j \m R ) IV 



cosh r 
1 



C N {rf 



(24) 



(25) 
(26) 
(27) 



for (EH), (22) and ([24]). 

On the other hand, the density matrices for the indi- 
vidual subsystems ([15]). (gDJ) , (ST) are 



— Y 



V tanh 2 ™ r 



Cjv(r) 2 ^ cosh 2 r 

v y n— 



sinh r 



\n)(n\ (28) 



pARR 


= I*X*I 


(15) 


p AR = 


Trivp ARR 


(16) 


p AR = 


Tr iP ARR 


(17) 


p«« = 


Tr M p ARR 


(18) 



and the density matrix for each individual subsystem 





= Tr iP AR = 


Tr IVP AR 


(19) 


P R 


= Tr Iv p RR 


= Tr M p AR 


(20) 


P R 


= Tr z p RR = 


Tr M p AR 


(21) 



p = 



p 

Where 



C N (r) 



'^tanh^r- 



cosh 2 r 

n— 



tanh 2JV r 
cosh r 



\N){N\ 



1 



cosh 2 r , 



\n)(n\ 



C N {rf 



D N{r) = y ^Jll = i _ (tannr) 2 ^ 1 ' 
^— 1 cosh r 



(29) 



(^(r)|0)(0|+^(r)|l)(l|) (30) 



(31) 
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N-l 



£>jv(r)=E( n+1 ) 



tanh 2 ™ r 



n=0 



cosh 2 ' 



= 1- 1 



N 



cosh 2 ', 



tanh 2 r 



(32) 

Notice that D^fr) + D^(r) = Cjy(r) and consequently 
all the density matrix traces are 1 as it must be. As all 
the probability is within the modes that we are consider- 
ing, all the possible decoherence is confined to the finite 
occupation number Hilbert space we are studying. 

As an effect of the imposition of the finite dimension 
PA —> |0)(0| asa->oo for any finite N, but it tends to 
i(|0><0j + |1)(1|) when A -> oo for all a. It is important 
to notice that both limits do not commute. The limit 
A — > oo should be taken first in order to recover the 
standard scalar field result. 
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. 71=0 



tanh r 
cosh 2 r 
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On (On + - — ; — tanh r 
coshr 

n + 1 



(|0n) 



)7Z + 1 
+ 5- |ln)(ln 
COsh r 
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+ ■ 



tanh 2JV r 
cosh 2 r 



| OA) (OA | 



C N {rf 



(35) 



V tanh™ +m r 



E 

m=0 



cosh 2 r 



\nm 



){mn\ + 



n=0 
m=0 



tanh" +T "r 
cosh 4 r 



x \/n + l\Jm + 1 1 71 + 1 m) (to + 1 7l| 



(36) 



IV. ANALYSIS OF CORRELATIONS 



In the following subsections we shall compute the neg- 
ativity of each bipartition of the system. 



In this section we will analyse the correlations tradeoff 
among all the possible bipartitions of the system. We will 
account for the entanglement by means of the negativ- 
ity, and we will study the total correlations by means of 
the mutual information, which accounts for both classical 
and quantum. 



1. Bipartition Alice-Rob 

Apart from the diagonal elements corresponding to 
|00)(00| and |1A)(1A| (which form two 1 x 1 blocks 
themselves), the partial transpose of the density ma- 
trix p AR (IMl) has a 2 x 2 block structure in the basis 
{lOn+lhlln)}^ 1 



A. Quantum correlations 



tanh 2 " r 



C(r) 2 cosh 2 



/tanh 2 , ^+1 
cosh r 
\Jn + 1 n 

V coshr sinh 2 r 



(37) 



Negativity is an entanglement monotone defined as the 
sum of the negative eigenvalues of the partial transpose 
density matrix for the system, which is defined as the 
transpose of only one of the subsystem q-dits in the bi- 
partite density matrix. If cr, are the eigenvalues of p p ^ B 
then 



a; 



AB 



E °* 



(33) 



ai<0 



To compute the negativity we need the partial trans- 
pose of the bipartite density matrices (|2"2"j) . (1231) and ([2"4"|) . 
which we will notate as r] AR , rj AR and r] RR respectively. 



Hence, the eigenvalues of are 

( n 



A 



tanh 2 ™ r 
2CV(r) 2 cosh 2 r 



V sinh 2 



tanh r 




tanh 



CAr(r) 2 cosh r 



vv+i 



cosh r 
A (tanh r) 



N-l 



n=0 
2N-2 



Cn (r) 2 cosh r 



(38) 



Where the notation \n= a means that n takes all the in- 
teger values from a\ to a at. 

Therefore the negativity for this bipartition is 



„AR 



tanh 2n r 



E 

. n=0 



cosh r 



| On) (Or 



Vn + 1 
coshr 



|0n+l)(lrc| 



)Tl -\- 1 
+ — — |ln + l)(In+l| 



cosh 



tanh 2JV r 
cosh 2 r 



|0A)(0A| 



C N {rf 



(34) 



N-l , , 2n 

M AR = Y '" 



■f Q 2CV(r) 2 cosh 2 r 



/ n 



V sinh 2 ', 



tanh r 



. sinh 



tanh 



cosh 



(39) 



Figure [5] shows the behaviour of negativity for all val- 
ues of A, which is clearly similar for all cases no matter 
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how many dimensions we are allowing for each mode. De- 
spite this, negativity AR is slightly sensitive to dimension 
variations. This is a difference with the fermionic cases, 
in which varying the dimension does not affect A/'f^ lionic 
at all in [23 |. 

M AR is shown in figure [3] as a function of r for differ- 
ent values of N, comparing them with the case N = 1. 
A very interesting result emerges here, for any pair of 
values for the maximum occupation number Ni < N2 
both negativity curves cross in a point a = a^Nx, N 2 ). 
This means that for any finite value of the Hilbert space 
dimension there is a region a < a c (N%, N 2 ) (low accelera- 
tions) in which entanglement is more degraded for higher 
dimension, and another region a > a c (Ni,N 2 ) (high ac- 
celerations) in which entanglement is more degraded for 
lower dimension. 

This disagrees the naive argument that higher dimen- 
sion would lead to higher Unruh decoherence which is 
not necessarily true. Figure 2] shows the behaviour of 
r c (l,N) as N grows. The crossing point with the neg- 
ativity curve for N — 1 grows as we consider larger N 
curves. r c (Ni, N2) is related with a c (Ni,N2) by means 
of the relationship 

For the limit N 2 — > 00, a c (Ni, N 2 ) — > 00, which means 
that infinite dimension negativity is below all the finite 
dimensional curves. 




r=arctanh(exp(-jt k c/a)) 



Figure 2: (Colour online) Bundle of curves displaying the 
negativity for the bipartition AR for all values of TV as a 
function of the acceleration 
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Figure 3: (Colour online): Same as fig. [5] for selected choices 
of TV showing the existence of crossing points. To the right 
(left) of these points entanglement degrades less (more) for 
lesser N. Solid blue line N — 1, dashed red line N = 2, 
dotted purple line N = 15, black dash-dotted line N = 00. 
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Figure 4: (Colour online) Values of r in which negativity 
curves for occupation number N cross the negativity curve 
for N = 1. In the region above the curve entanglement de- 
grades faster for N — 1 than for N > 1. 



2. Bipartition Alice-AntiRob 



{lOnKlln + l}}^ 1 



Excepting the diagonal elements corresponding to 
1 10) <10| and |0iV)((W| (which form two 1 x 1 blocks 
themselves), the partial transpose of the density ma- 
trix p AR (|55|) has a 2 x 2 block structure in the basis 



tanh 2 ™ r 
C(r) 2 cosh 2 



tanhr 
\ cosh r 



Vn+l 



tanln 
cosh? 
tanh 2 ?' 



cosh 2 r 



Vn + T \ 



(n + 2) 



(40) 
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Hence, the eigenvalues of (|35l) are 



tanh 2 " r 
2C(r) 2 cosh 2 



1 + (n + 2) 



tanh 2 r 
cosh 2 r 



± 



1 + (n + 2) 



tanh 2 
cosh 2 ; 



/ \ 4 tanh 2 r 



cosh 2 : 



N-l 



71=0 



X 



N 



C(r) 2 cosh r ' 



A 



jv+i = 



tanh 2W r 
C(r) 2 cosh 2 r 



(41) 



Therefore, the negativity for this bipartition is always 0, 
independently of the value of the acceleration parameter 
and the occupation number bound N. 

From this results can be concluded that limiting the 
dimension has no effect in the creation or not of quantum 
correlations between Alice and AntiRob. As far as the 
field is bosonic, no entanglement is created in the CCA 
bipartitions of the system no matter how we limit the 
dimension of the Hilbert space. 



3. Bipartition Rob-AntiRob 



The partial transpose of the density matrix p RR (|36|) 
has a block structure. Namely, it is formed by 2N + 1 
blocks whose dimension varies. In the following we will 
detailedly analyse the blocks. 

1. First of all, we have N + 1 blocks {M fl }^} which 
are endomorphisms that act in the subspace (of di- 
mension D) expanded by the basis Br) = {|mn)} 
in which to + n = D — 1 < N . 

2. Then we have N more blocks {M' D }^ =1 that 
act in the subspace (of dimension D) ex- 
panded by the basis B' D = {\m'n')} in which 
to' + n' = 2N - D + 1 > N. Notice that not all 
the possible to' and n' are allowed due to the limi- 
tation to the occupation number m',n' < N. 

As an example which will perfectly clarify this 
construction, if N = 4 there will be 9 blocks, 
Ml, M 2 , M 3 , M 4 , M 5 , Mi, Mg, M' 2 ,M' X each one is an en- 
domorphism which acts in the subspace expanded by the 
bases 



B l 


= {|oo)} 






B 2 


= {|oi>, 


|io>} 




B, 


= {102), 


120), 


11)} 


Bi 


= {|03>, 


130), 


12),|21)} 


B 5 


= {|04>, 


140), 


13),|31),|22)} 


B' 4 


= {|14>, 


|41), 


23),|32)} 


B' 3 


= {|24>, 


142), 


33)} 


B' 2 


= {|34>, 


|43)} 






= {|44>} 







(42) 



respectively. 

In this fashion, the whole matrix is an endomorphism 
within the subspace R=0^ = H L 1 >%©0j! =1 S'p being Si the 
subspace (of dimension D — i) expanded by the basis Bi 
and S'j the subspace (of dimension D = j) expanded by 
the basis B'j . 

The blocks Mi,...,Mjv+i and M[, 
form the matrix J36[) have the following form 



M' N which 



M 



D 



(o 


a i 










\ 


ai 





«2 













«2 





03 













«3 





04 ■ 


















•. ••• 
















a D -i 


\o 











. 


. . a D -i a D J 














\ 


bi 





b 2 







.... 





b 2 





63 




.... 








63 





b 4 ■ 


.... 















. ••• 
















■ b D -i 


V° 











. 


■ b D -i b D J 



(43) 



(44) 



The matrix elements a n and b n are defined as follows 
(tanhr) 13 - 1 



02;+! = 



C(r) 2 cosh r 

VD-iVi- {t&nhr)D ~ 2 



921+1 



C(r) 2 cosh r 
(tanhr) 2jv - g+1 
C(r) 2 cosh 2 r 



b 2 i = + 1-WI + N -D + l 



(tanh i 



\2N-D 



C(r) 2 cosh r 

(45) 



Notice that the elements arc completely different when 
the value of the label n is odd or even. 

As the whole matrix is the direct sum of the blocks 



'N+l 



N 



RR 



(46) 



\D=l 



\D=l 



the eigenvalues and, specifically, the negative eigenvalues 
of r] RR would be the negative eigenvalues of all the blocks 
Mr) and M' D gathered togheter, which can be easily com- 
puted numerically. Figure [5] shows the behaviour of AT RR 
with r and for different values of N. 

We can now compare the finite N bosonic case with 
their same dimension analog for fermions. Namely, a 



8 




r=arctanh(exp(-7l k Q c/a)) 



Figure 5: (Colour online) Negativity RR for different values 
of TV, showing the upper bound reached when a —¥ oo. Nega- 
tivity diverges in the black hole limit only for TV — » oo. 

Grassman scalar field (spinless fcrmion) has the same 
Hilbert space dimension as the scalar case with TV = 1, 
the only difference is the anticommutation of the field 
operators instead of the commutation which applies for 
bosons. On the other hand, scalars limited to TV = 4 
and TV = 2 can be considered as two different analogs to 
the Dirac field as the former has the same Hilbert space 
dimension as Dirac modes and the latter would share the 
same possible maximum occupation number. 

This comparison can be seen in figures [6j [71 We see 
that the comportment is similar (monotonic growth from 
zero to a finite limit for a — > oo) but the functional de- 
pendence is still very different in both cases. Specifically, 
as a increases the bosonic cases grow a higher entangle- 
ment between the modes of the field on both sides of the 
horizon than the same dimension fermionic analogs. The 
ability to create correlations between modes in / and IV 
is related with the dimension of the Hilbert space but not 
only. Statistics plays also an important role. 

This clearly shows another important difference be- 
tween fermionic and bosonic fields. Pauli exclusion prin- 
ciple prevents the total degradation of fermionic entan- 
glement in CCA bipartitions, whereas, conversely, im- 
pedes entanglement creation between RR. 

B. Mutual information 



0.6h 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

r=arctan(exp(-ji k c/a)) 



Figure 6: (Colour online) Negativity of the bipartition RR 
for fermion fields. Namely Grassman scalar (red dashed line) 
and Dirac (blue solid line). Negativity upper bound is greater 
for the Dirac case as dim('HDirac) > dim(%Grassman.) Notice 
that here r = ta,n(e~ 7Tk °^ a ) instead of the hyperbolic tangent 
and therefore r — > 7r/4 =>■ a oo. 



2r 

1.8- Scalar N=4 

1.6- / 




Scalar N=2 



Scalar N=1 



r=arctanh(exp(-7l k Q c/a)) 

Figure 7: (Colour online) Same as fig. [6]but for bounded occu- 
pation number scalar fields. TV = 1 (red dashed line) is dimen- 
sionally analogous to the Grassman scalar case. TV = 4 (black 
dash-dotted line) is dimensionally analogous to the Dirac case. 
TV = 2 (blue solid line) is analogous to the Dirac field in max- 
imum occupation number. Notice that more entanglement is 
grown for greater TV. 



Mutual information accounts for correlations (both 
quantum and classical) between two different partitions 
of a system. It is defined as 

Iab =S a + S b - S ab (47) 

where Sa, Sb and Sab are respectively the Von Neu- 
mann entropies for the individual subsystems A and B 



and for the joint system AB. 

To compute the mutual information for each biparti- 
tion we will need the eigenvalues of the corresponding 
density matrices. We shall go through all the process 
detailedly in the lines below. 
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1. Bipartition Alice-Rob 

Excepting the element | OAT) (OAT | (which forms a 1 x 1 
block itself) the density matrix for the system Alice- 
Rob ([22jl consists on JV 2 x 2 blocks in the basis 
{|0n) , |1 n + l)}n=o which have the form 



tanh 2 " : 



1 



Vn+ 1 

cosh r 

C N {r) 2 cosh 2 r \ V n + 1 n + 1 



coshr cosh z r 
Hence, the eigenvalues of (|22[) are 

n + 1 



An 



A 



tanh 2 ™ r 
CV(r) 2 cosh 2 r 
tanh 27V r 



1 



cosh 2 r 



N 



CV(r) 2 cosh 2 r 
Along with N identically zero eigenvalues. 



N-l 



n=0 



(48) 



(49) 



It is easy to see that rank(X) = rank(Y") = 1. This 
means that all the eigenvalues of (j2"4"]l are zero except for 
two of them, which we can readily compute 



\RR 



C(r) 2 



\RR 

Ay 



C{r) 2 



(54) 



where D N (r) and D]^(r) are given by (|3Tj) and (|32 



4- Von Neumann entropies for each subsystem and mutual 
information 



To compute the Von Neumann entropies we need the 
eigenvalues of every bipartition and the individual den- 
sity matrices. The eigenvalues of p AR , p AR , p RR are 
respectively (|39"]l. (fSTjl and ([54]) . 

The eigenvalues of the individual systems density ma- 
trices can be directly read from (|2"5f , (|2"9")) and (f5U|) since 
p fl and p A have diagonal forms in the Fock basis. The 
Von Neumann entropy for a partition B of the system is 
S = — Tr(plog 2 p). Therefore their entropies are 



2. Bipartition Alice-AntiRob 

Except from the diagonal element corresponding to 
1 00) (00 1 (which forms one lxl block itself) the den- 
sity matrix for the system Alice-AntiRob (|23|) consists 
on V - 1 2 x 2 blocks in the basis {|0n) , \ln - 1)}„ =1 
which have the form 



tanh 2 ™ r 
Cn{t) 2 cosh 2 r 

Therefore the eigenvalues of 



1 



sinh r 

n 



(50) 



arc 



tanh 2 ™ r 



Cn{t) 2 cosh r 



sinh r 



N 



n=0 



(51) 



along with V identically zero eigenvalues. 



3. Bipartition Rob-AntiRob 

The density matrix for Rob-AntiRob consists in 
the direct sum of two blocks 



p RR = X®Y 



(52) 



of dimensions dim(V) = N+1, dim(F) = V. The matrix 
elements of X and Y are 



Xij — 



(tanh r) 



i+j-2 



C'n (r) 2 cosh r 



Y, 



ij — V l v J 



(tanh r) 



i+j-2 



Cn (r) 2 cosh r 



(53) 

in the bases {\nn)} n=Q and {\n + ln)} n=Q respectively. 



N , ,2n 

tanh 



Cn{t) 2 cosh r 



x lo 



2 

N-l 



tanh 2 " r 



CW (r) 2 cosh r 



sinh r 
n 



S n 



E 



tanh 2 ™ r 



q Cn(t) 2 cosh 



sinh r 
n + 1 



x log 2 



tanh 2 ™ r 



Cn (r) 2 cosh r 



Si 



tanh 2JV r 
CAr(r) 2 cosh 2 r 

21og 2 [Cjv(r)] 



log 2 



cosh 2 r 

n + 1 

H o— 

cosh r 

tanh 2W r 
Cn(t) 2 cosh 2 r 



Cjv(r) 



J] ^(r)log a [D N (r)] 
(55) 



We obtain a universal result which relates the entropies 
of the different bipartitions of the system, 

Sar = Sr, Sar = Sr> Srr — $a (56) 
These results can be summarised in the expression 



Sij = Sr 



(57) 



where /, J and K labels represent different subsystems. 
Whichever values I ^ J ^ K will satisfy the identity. 
This is also true for standard scalar fields, Grassman 
scalar fields (spinless fermions) and Dirac fields. These 
relationships are completely universal, being independent 
of statistics and dimension, and reflect a fundamental as- 
pect of Unruh decoherence in terms of the entropy of the 
partial systems, namely, the way in which the entropy of 
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the bipartitions behaves as acceleration increases is not 
independent from the way the individual entropies do. 

The mutual information for all the possible bipartitions 
of the system will be 

Iar = Sa + Sr — Sar = Sa + Sr — S r 
Iar. — Sa + Sr — Sar — Sa + Sr — Sr, 
Irr — Sr + Sr — Srr — Sa + Sr + Sr 

The first notable difference from the standard bosonic 
field is that we do not obtain here the conservation law 
of the mutual information for the system Alice- Rob and 
Alice- AntiRob 

Iar + I A R = 2S A (58) 

And that for any finite N, Sa goes to zero when a —> oo. 
Only in the limit of iV — > oo the conservation law is 
restored. 

Figure [5] shows the mutual information tradeoff of the 
systems AR and AR from N = 1 to N = 10 4 , along 
with the limit N — > oo in which the conservation law is 
fulfilled for all values of a. 



0.2- \ \ \ \ \ 

s N V v x 

I 1 " - ~ - - - — = = -V - - - -- 1 

2 4 6 8 10 

r=arctanh(exp(-7l k Q c/a)) 

Figure 9: (Colour online) Violation of mutual information 
conservation law for the systems AR and AR. The conser- 
vation law is fulfilled until a reaches a critical value which 
is logarithmically displaced to the right as N increases. The 
conservation law is completely restored when TV — > oo 



\ Alice-Rob 
















X S \ N=10 2 \ 




N=1o\ 




*-A N=1 °\ 

» N \ \ 
* \ \ 

N=1 \ \ 
1 Xy. N 

/ Alice-AntiRob 

1 1 1 T 









2 3 4 5 

r=arctanh(exp(-ji k c/a)) 




Figure 8: (Colour online) Mutual information for the sys- 
tems Alice-Rob (blue continuous lines) and Alice-AntiRob 
(red dashed lines) as the acceleration parameter varies. Sev- 
eral values of A^ are plotted along with the N = oo case. A 
conservation law is satisfied until acceleration reaches a criti- 
cal value which is displaced to the right logarithmically as N 
increases. 

As it can be seen in figure |H] the largest deviation from 
the conservation law is obtained for N = 1. As it is 
shown in the figure, for a given N the conservation law 
is fulfilled until the acceleration reaches a critical value 
a = ai, then correlations go rapidly to zero. This critical 
value increases logarithmically with N. 

We showed above that quantum correlations between 
AR are quickly lost as a increases for all N and no entan- 
glement is created between AR. This means that for the 



Figure 10: (Colour online) Behaviour of the mutual informa- 
tion for the system Rob-AntiRob as acceleration varies for 
different values of Af showing the upper bound reached when 
a — s> oo. Mutual information diverges in the black hole limit 
only for N — s> oo. 



high acceleration regime (where all quantum correlations 
vanish) classical correlations dominate mutual informa- 
tion. Therefore, what we learn from mutual information 
in this regime is the behaviour of purely classical cor- 
relations which are usually very difficult to be studied 
separately from the quantum ones. 

For the scalar field we have seen that, conversely to 
fermionic fields [13, [HI , limiting the dimension produces 
boundary effects which make classical correlations go to 
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zero. Here, conservation of these correlations for all val- 
ues of a requires infinite dimension. In other words, finite 
dimensions schemes kill classical correlations as Rob ac- 
celerates. 

One would expect something similar for fermions since 
their states are naturally of finite dimension. Hence, sim- 
ilar 'border effects' in classical correlations should appear 
in the same fashion as for bosons. Despite this fact, mu- 
tual information for fermions does not vanish. The ex- 
planation for this difference between bosons and fermions 
comes from fermionic quantum correlations. As shown in 
there is a conservation law for fermionic quantum 
entanglement 

KfAR , KfAR _ }_ (rq\ 

J v fermions * J v fermions ^ \°^) 

Since classical correlations for the finite dimensional case 
go eventually to zero, quantum correlations rule mutual 
information behavior for fermions. Therefore, it can be 
concluded that the conservation law for the mutual in- 
formation for fermions must be strongly related with the 
conservation of the fermionic entanglement which has his 
origin in statistics. 

We can conclude then that the origin of the univer- 
sal mutual information conservation law is different for 
fermions and bosons. On one hand, for bosons, it ap- 
pears as a classical correlations conservation law. On 
the other hand, for fermions, this law reflects a quantum 
correlations conservation. This can also explain why mu- 
tual information behaves so similarly to negativity for 
fermions as it was obtained in [23[. 

To finish this work and complete the analysis of mutual 
information let us show in figure [10] how the behaviour 
of I R ft changes as N is incresaed and how the diver- 
gent limit is obtained when N — ¥ oo. The results about 
mutual information here are coherent with the thorough 
analysis of the correlations for the RR bipartition per- 
formed above when we analysed negativity. 

V. CONCLUSIONS 

In this paper we answer the question of the actual im- 
pact of Fock space dimensionality on Unruh decoherence 
phenomena. To do so, we have studied the dimensional 
dependence of scalar field correlations when one observer 
is non-inertial. With this end in sight we have built a 
scalar field entangled state in which we have imposed a 
maximum occupation number N for Rindler modes. 

We have shown that the comportment of the entan- 
glement for AR and AR is only slightly influenced by N 
so in other words, the qualitative behaviour (quick loss 
of entanglement for AR as shown in figure [2] and no en- 
tanglement creation for the system AR) is the same for 
finite and infinite N. This again points to the argument 
in [22|, that it is statistics and not dimensionality that 
conditions the comportment of correlations in the pres- 
ence of horizons. 



However, we have shown that AR entanglement is sen- 
sitive to variations of N. This opposes what we found 
for fermions, whose correlations are completely insensi- 
tive to Hilbert space dimension variations (for example 
going from Grassman scalars to Dirac fields). 

In previous works we found a universal behaviour for 
the fermionic entanglement of the bipartition AR. Specif- 
ically the functional form for the negativity was exactly 
the same; independent of the maximally entangled state 
selected, the spin of the field, and the number of modes 
considered going beyond SMA. In all the cases Unruh de- 
coherence degrades fermionic entanglement exactly the 
same way. Here we see that for bosons this universality 
principle does exist but it is not as strong due to the 
sensitivity of AR to dimension changes. Extending this 
formalism to the electromagnetic field will reveal more 
information about this, thus it is expected to appear else- 
where. 

We have also seen that lesser N does not necessarily 
imply faster entanglement degradation. Instead, we have 
shown that for two different finite values of N, namely 
N\ < N2 there is a region a < a c in which entanglement 
is more degraded for N2 and another region a > a c in 
which entanglement is more degraded for N\. In other 
words, for high accelerations, higher dimension means 
less entanglement degradation by Unruh effect. This re- 
sult clashes again with the extended idea that lesser di- 
mension would protect correlations better than higher 
dimension, one misconception that after all this research 
should be banished from the explanation of these phe- 
nomena. 

We have also showed that, since a c shifts to the right 
as N is increased, in the limit N — > 00, a c — > 00, so 
that entanglement is more degraded for the infinite di- 
mensional case than for any finite N whatever the value 
of the acceleration. 

It is remarkable that there is no entanglement creation 
in the CCA bipartitions even for finite dimension; no en- 
tanglement is created for the bipartition AR whatever 
the dimension limit N. This reflects again that the dif- 
ferences between fermions and bosons have nothing to 
do with the finite dimensionality, but with the different 
statistics. 

Concerning mutual information, we have shown that 
the conservation law found for scalar fields [23[ for the 
systems AR and AR is violated for finite values of N. 
We have obtained that for a finite N the conservation 
law is fulfilled until acceleration reaches a critical value, 
in which correlations quickly drop. This critical value 
grows logarithmically with N, which means that the con- 
servation law is satisfied for all a when N — > 00. There- 
fore, the violation of this conservation law can be associ- 
ated with the boundary effects of imposing a dimensional 
limit. It is important to observe that in the bosonic case, 
mutual information is mainly accounting for classical cor- 
relations. This is because quantum correlations among 
these systems are quickly lost as a increases. 

Surprisingly, fermions, which have necessarily a lim- 
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ited dimension Hilbert space for each mode, lack these 
boundary effects. This big difference between fermions 
and bosons is related with the statistical conservation of 
quantum correlations in the black hole limit for fermions. 
Conversely to bosons, the high acceleration regime for 
fermions does not mean that mutual information is re- 
flecting any classical correlations. Instead, mutual infor- 
mation is mainly reflecting fermionic quantum correla- 
tions. These correlations satisfy themselves a conserva- 
tion law (1591) which is 'inherited' by mutual information. 
This also explains the similitude between mutual infor- 
mation and negativity behaviour for fermions found in 
23]. 

The conclusion here is that the universal conservation 
law for mutual information is found for both fermions and 
bosons, however the nature of this conservation is com- 
pletely different. With reference to bosons it is due to 
classical correlations conservation, whereas for fermions 
it is due to quantum correlations conservation. This il- 
lustrates yet again that statistics is such a paramount 
feature in order to explain how correlations behave in 
the presence of horizons. 

The dimension of the Fock space has the largest impact 
in the behaviour of correlations between regions I and 
IV separated by the horizon. Comparing the limited di- 
mension scalar fields with their fermionic analogs we have 
found that the behaviour of correlations RR is somewhat 
similar for bosons and fermions and mainly ruled by the 
dimensionality of the Hilbert space. For both fermionic 
and bosonic fields, entanglement is always created be- 
tween RR reaching a maximum value when a — > oo for 



any finite dimension. 

The scalar cases of finite N present, however, impor- 
tant differences compared with their fermionic Fock space 
dimensional analogs (namely N = 1 is analogous to the 
Grassman scalar case and TV = 2,4 to the Dirac field 
case). In the black hole limit, scalar states entanglement 
created between modes in both sides of the horizon is 
greater than in the fermionic case. 

The scalar cases of finite N however present important 
differences compared with their fermionic Fock space di- 
mensional analogs (namely N = 1 is analogous to the 
Grassman scalar case and N — 2,4 to the Dirac field 
case). In the black hole limit the scalar states entangle- 
ment created between modes in both sides of the horizon 
is greater than in the fermionic case. 

This implies that, for correlations between modes of 
Rob and AntiRob, the effect of statistics in the black 
hole limit is the opposite to what happened for the CCA 
bipartitions. Here bosonic correlations reach higher val- 
ues than their corresponding fermionic analogs. 
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